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Abstract: We compute direct current (dc) thermoelectric transport coefficients in strongly coupled 
quantum field theories without long lived quasiparticles, at finite temperature and charge 
density, and disordered on long wavelengths compared to the length scale of local thermaliza- 
tion. Many previous transport computations in strongly coupled systems are interpretable 
hydrodynamically, despite formally going beyond the hydrodynamic regime. This includes 
momentum relaxation times previously derived by the memory matrix formalism, and non- 
perturbative holographic results; in the latter case, this is subject to some important sub¬ 
tleties. Our formalism may extend some memory matrix computations to higher orders in 
the perturbative disorder strength, as well as give valuable insight into non-perturbative 
regimes. Strongly coupled metals with quantum critical contributions to transport gener- 
ically transition between coherent and incoherent metals as disorder strength is increased 
at fixed temperature, analogous to mean field holographic treatments of disorder. From a 
condensed matter perspective, our theory generalizes the resistor network approximation, 
and associated variational techniques, to strongly interacting systems where momentum is 
long lived. 
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Introduction 

One of the most exotic and mysterious systems in condensed matter physics is the strange metal, non- 
Fermi liquid phase of the high Tc superconductors [1, 2]. The transport data in these materials - including, 
most famously, the linear in temperature dc electrical resistivity -- defies clear explanation by a theory of 
long lived quasiparticles [3]. Alternatively, the effectively relativistic plasma in graphene may provide an 
experimental realization of a strongly interacting quantum fluid [4]. Finally, recent advances in ultracold 
atomic gases have paved the way to realizing strongly interacting fluids [5]. In all of the above systems, the 
absence of long lived quasiparticles on experimentally appropriate time scales (e.g., in the computation of 
dc transport coefficients) poses a challenge for traditional, quasiparticle-based approaches to condensed 
matter physics. 

From a theoretical perspective, a generic strongly interacting quantum field theory (QFT) in more 
than one spatial dimension has only a few quantities (energy, charge and momentum) that are long lived, 
and so hydrodynamics may be a sensible description of the low energy physics at finite temperature 
and density of all of the above systems. Though hydrodynamics is an old theory [6, 7], its implications 
for transport have only been understood comparatively recently [8], in weakly disordered systems near 
quantum criticality, in external magnetic fields [8, 9, 10], and in some simple examples of disordered, 
non-relativistic electron fluids [11]. This is because “textbook” hydrodynamics is utterly inappropriate 
for most metals, where the electron-impurity scattering length is short compared to the electron-electron 
scattering length. Momentum and energy rapidly decay, and the only hydrodynamic variable is the charge 
density. Note that in contrast to this canonical lore, [12] proposed that observing viscous hydrodynamics 
in some metals may be possible. 
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In most ways, hydrodynamics is a far simpler theory to understand (and perform computations in) than 
quasiparticle based approaches, such as kinetic theory. The difficulty in studying these systems theoreti¬ 
cally lies in the fact that hydrodynamics does not completely solve the transport problem: the coefficients 
in the hydrodynamic equations must be related to Green’s functions in a microscopic model. Nonetheless, 
if hydrodynamics is valid, it does provide universal constraints on transport, and a transparent physical 
picture to interpret the results. There are two tractable approaches that can compute the requisite mi¬ 
croscopic Green’s functions, without reference to quasiparticles. The first is methods from (perturbative) 
QFT, combined with the memory matrix approach [13, 14, 15], which has recently been used in many mi¬ 
croscopic models of strange metals, reasonable for describing cuprate strange metals [16, 17, 18, 19]. These 
approaches rely on properly resumming certain families of Feynman diagrams to all orders. The second 
approach is holography [20, 21, 22], which reduces the computation of Green’s functions to solving classical 
differential equations in a black hole background. This can be done numerically [23, 24, 25, 26, 27, 28, 29], 
though in some cases analytic insight can be obtained [30, 26, 31, 32, 28, 33, 34, 29, 35], sometimes by 
employing the memory matrix method [8, 36, 37]. 

Surprisingly, many of the above transport theories from recent years completely match hydrodynamic 
predictions, at least superficially, despite being formally beyond the regime of validity of hydrodynamics. 
We take this as an indication that a thorough understanding of hydrodynamic implications for transport 
in disordered theories is worthwhile, though we will also carefully describe the regime of validity of 
the approach. In addition, while almost every citation above aims to address the strange metal phase 
[1, 2, 3], the “hydrodynamic insight” gained from these methods may be applicable to a much broader 
set of experimentally realized interacting quantum systems. 

1.1 I Motivation: Incoherent Metals and Holography 

Let us begin with the main quantitative motivation for the present paper, which is the physical interpre¬ 
tation of a large body of recent holographic work on transport in QFTs without translational symmetry. 

[8] proposed a simple hydrodynamic framework for dc transport which has been quite predictive of 
both holographic and memory function results in subsequent works, at weak disorder. As we previously 
mentioned, this framework has been surprisingly good at describing many holographic models that treat 
disorder at a mean field level. A natural conjecture is that disordered hydrodynamic dc transport can 
describe holographic systems with explicitly broken translational symmetry, and so it is worthwhile to 
fully flesh out the disordered hydrodynamic formalism. 

We begin with a generic hydrodynamic framework for zero frequency transport calculations in Section 
2. Our emphasis is on a clear presentation of the assumptions and regime of validity of a hydrodynamic 
description of transport. 

We exactly solve the transport problem in Section 3, at leading order in the strength of disorder, in 
the limit where translational symmetry is only broken weakly. These systems describe coherent metals, 
in the language of [38] - henceforth we will also adopt this terminology. We show that our resulting 
computations exactly equal the results found by the memory function formalism, under the assumption 
that the momentum is long-lived, justifying that our approach is sensible, as well as providing a physically 
transparent derivation of memory function based formulas for conductivities (at least, in the mutual regime 
of validity of the two methods). 

We further show in Section 4 that the hydrodynamic framework can be used to interpret exact, non- 
perturbative analytic results for dc transport found using holography. This is subject to the important 
subtlety that the transport coefficients are computed in terms of a new emergent horizon fluid with 
distinct, emergent (but somewhat sensible) equations of state. 

We then proceed to study hydrodynamic transport in non-perturbatively disordered QFTs. Though 
not amenable to analytic techniques, we develop a combination of rigorous variational approaches and 
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heuristic approximations, outlined in Section 5, to calculate dc transport in this regime. One might expect 
that transport becomes dominated by dissipative hydrodynamics, as momentum may become a “bad” 
conserved quantity; such a state is an incoherent metal, in the language of [38]. We find further evidence 
for this qualitative picture. 

To date, all models of incoherent metals are holographic massive gravity models [39, 40, 41, 42, 43, 
44], or similarly inspired holographic approaches [45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55] which we 
will henceforth lump together under the label of mean-field disorder. These models break translational 
symmetry phenomenologically, but not explicitly.^ These models always predict dc transport which, 
at all disorder strengths, can be interpreted in terms of the hydrodynamic results of [8], or the slight 
generalization of [56]. A simple example of this is the exact formula for dc electrical conductivity in an 
isotropic system [41]: 


a = aQ + 


Q^t 
M ’ 


( 1 ) 


with (Tq a transport coefficient independent of disorder strength, Q a charge density, and Af an analogue 
of the mass density. The parameter r is analogous to a momentum relaxation time, and related to the 
phenomenological graviton mass. This formula was already known from quantum critical hydrodynamics 
[8], using computations valid as r —>■ oo. Indeed, the latter term is nothing more than the Drude formula, 
valid in a system without quasiparticles, and the former is a quantum effect that can be important close 
to a particle-hole symmetric point [57]. Mean field models always predict that (1) holds even as r —?■ 0, or 
in the non-perturbative, strong disorder regime. In this limit r cannot be interpreted as the momentum 
relaxation time directly, but importantly, cr stays larger than ctq, which is the conductivity when Q = 0 
(an uncharged theory).^ And while mean field models do agree with approaches that explicitly break 
translational symmetry weakly [31, 32, 33], this is simply a consequence of the perturbative equivalence 
between holographic and memory function computations of transport, proven in many cases in [58]. 

One might suspect that the fact that (1) holds as r —)■ 0 is a sign that mean field physics is a 
poor description of a strongly disordered QFT, even in holography. For example, it is well known that 
mean field descriptions of disorder can completely fail to capture even basic thermodynamics of strongly 
disordered spin models in classical statistical physics - instead, the emergent phases are spin glasses and 
must be treated using much more delicate technologies [59]. 

Our work in Section 5 demonstrates that our hydrodynamic framework gives an independent frame¬ 
work in which the qualitative picture of dc transport given in (1) is correct at all disorder strengths until 
the hydrodynamic description fails. As an important example, we argue that for an isotropic quantum 
critical system where viscous transport may be neglected, 


0? 


crQi(u) -h (Tq2(u)^ < a 


< 


Qn 

0-Q3(w) +<TQ4(n)^, 


( 2 ) 


with a the dc electrical conductivity, Qo the spatial average of the charge density Q, u the typical size 
of fluctuations in Q about this average, and crQi, 2 , 3,4 are related to the “quantum critical” diffusive con¬ 
ductivity, cjq.^ As m —)• 0, they are all proportional to the constant ctq associated with the translationally 
invariant QFT. At stronger u, <Tqi^ 2 , 3,4 may be complicated (spatially-averaged) correlations between Q 
and (TqI see (91), (110) and (126). 

We have written (2) in the manner we did to emphasize asymptotic behavior. When u <C Qo) ^ ~ 
if Qo 7 ^ 0, and this is a direct consequence of the fact that currents can only decay due to the very slow 
relaxation of momentum, as in a translationally invariant (momentum conserving) system at finite Qo, 

^This is technically not quite right - there is one (set of) bulk scalar fields in these models which is of the form cj)i = kxi, 
but this choice maintains homogeneity in the sectors of the theory of interest. 

^See [54, 55] for recent updates on this particular holographic model. 

®In graphene, for example, Uq ~ e^/h, with e the charge of the electron. 
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Figure 1: A qualitative sketch of the coherent-incoherent transition realizable in our framework. 
a denotes the value of a transport coefficient, such as electrical conductivity, and u denotes the 
“strength of randomness”. The solid black line shows our perturbative analytic computation of 
a ~ as u —)• 0. The dashed red line is the qualitative prediction of mean field models that 

a saturates at a finite value at strong disorder in a theory with quantum critical transport; in 
particular, a > a*. The gray shaded region corresponds to the region of a allowed by variational 
bounds on a, in generic agreement with mean field models, tt ~ Qo is the scale of the crossover 
between a coherent and incoherent metal. 

Galilean invariance imposes a = oo.^ In contrast, when u 3> Qo, o' is sensitive to typical behavior of local 
cjq and is not parametrically larger. Remarkably, cjqi > 0 in any system where the local quantum critical 
conductivity never vanishes, so any such system is provably a conductor. The physical intuition behind 
this is that a current can always flow locally due to finite cjq, and so if these local currents can always 
flow, a global current flow can necessarily be established. The upper bound is simply a statement that 
(up to subtleties involving conservation laws) we can bound the power dissipated (and accordingly the 
conductance), with the average electric field fixed, by assuming that the electric field within the system 
is uniform. 

When u <C Qo) we can make contact with the memory matrix formalism and exactly compute a 
perturbatively, as we discuss in Section 3, and so in this regime we can do better than the bounds in 
(2). Indeed, in this regime, we can identify t/M. ~ aqlu^, and so our bounds justify (1) in our class 
of hydrodynamic models. Analogous phenomena may be responsible for the finite conductivity of mean 
field holographic models at all “disorder strengths”. 

A pictorial summary of (2) is shown in Figure I, and the main quantitative result of this paper is 
the justification for Figure 1 without any mean-field treatment of disorder, and the development of new 
techniques to address the strongly disordered regime. 

After we had completed this work, [60, 61] appeared, and has some overlap with ideas in Section 4. 

^We can maintain an electric current without adding any energy by simply shifting to a moving reference frame. 
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1.2 I Motivation: Beyond Resistor Lattices 

Though the main quantitative focus of this work is a set of computational tools to study hydrodynamic 
transport in relativistic fluids, such as in holography, we also emphasize that the framework we are 
developing (with suitable generalizations) is sensible for a description of transport in strongly interacting 
condensed matter systems, without any reference to holography. 

A common approximation made in condensed matter is what we will refer to below as the “resistor 
lattice” approximation, which in physical terms is the statement that the slow, hydrodynamic sector of 
the theory consists of only a conserved charge. One may then model the emergent hydrodynamics - a 
simple diffusion equation for charge - as a local resistor network: see e.g. [62, 63]. As mentioned before, 
this is sensible if electrons scatter more frequently from impurities than they do from each other. 

However, we will point out in Section 3.2 that this approximation fails in a clean hydrodynamic system: 
the necessary resistor lattice becomes nonlocal. This is not a surprise. What this paper clarifies is the 
technique that unifies the computation of transport in a weakly disordered (memory matrix) regime and 
a strongly disordered regime. In doing so, we generalize well-known variational techniques from resistor 
networks to account for convective transport. In very special cases [11] performed similar calculations, 
though did not elucidate the connections with the memory matrix formalism, or with resistor lattice 
technologies, which we generalize directly in the continuum in this paper. Such resistor lattice methods 
- commonly with an additional approximation called effective medium theory [64, 65] - have been used 
recently to study transport in a variety of experimentally realizable systems [66, 67, 68]. Our approach 
can generalize these computations to the regime when disorder is weak, and may result in interesting new 
experimental predictions. 

We emphasize that the calculations in [66, 67, 68] typically include non-relativistic effects such as 
Coulomb screening, or additionally approximate that electron-hole recombination is slow enough that 
both the electron and hole densities are hydrodynamic quantities. We will not make either assumption 
in this paper, but the general framework and many computational methods we develop almost certainly 
extend quite naturally to account for these effects. 



Steady-State Hydrodynamics 


We consider a strongly coupled QFT in d spatial dimensions at finite temperature and density, on a 
flat spacetime. It is necessary to generalize to curved spaces to connect with the results of Section 4, 
but every result in this paper generalizes in the obvious way (replacing partial derivatives with covariant 
derivatives, f d'^x —)• f d'^Xy^, etc.), and so we will not do so explicitly for ease of presentation. Without 
quasiparticles, the long time dynamics are that of charge, energy and momentum. In this section, we will 
work with relativistic notation, though the techniques work for non-relativistic theories as well. We focus 
on theories with a single conserved charge, but the techniques straightforwardly generalize to theories 
with multiple conserved charges. Note that we will work in units with h = 1. 

We deform the microscopic Hamiltonian H by an external chemical potential: 

H ^ H - j . (3) 

with a conserved electrical current, F = dA, and 

A = /2(x)dt, (4) 

so if Q(x) is the local charge density operator, 

/d'^x)2(x)Q(x). (5) 


6 



◄-► 



Figure 2: We employ a separation of 3 length scales in this paper, p., and the local fluid properties 
such as entropy density S, may vary substantially over the distance scale We require I <C C 
for a hydrodynamic description to be sensible. We will often put our fluids in a large but finite 
box of length L ^ as well. 


The chemical potential in the fluid is thus p. We also assume that the temperature is uniformly T, and 
that there is no fluid velocity, in our background state. This forms the basis of a consistent solution to 

hydrodynamic equations, driven by the coupling p to an external bath, as we will derive below. The 

steady-state hydrodynamic equations read (in relativistic notation) [8] 

^ (6a) 

dif = 0, (6b) 


where Greek indices denote spacetime indices and Latin indices denote spatial indices and is the 
energy-momentum current. We have implicitly taken expectation values over all operators in (6) and will 
do so for the remainder of the paper. Because we have sourced disorder in our fluid entirely through p{x), 
we do not need to couple any other dynamical sectors to the theory, though we will point out how this 
may be done perturbatively in Section 3.1, when additional scalars contribute to disorder. The coupling 
of the fluid to an external chemical potential means that both energy and momentum may be exchanged 
with the external bath. 

In order for hydrodynamics to be valid, it is necessary that p vary slowly in space, on a length 
scale ^ which is large compared to the (possibly position-dependent) mean free path of the fluid 1. In 
our strongly interacting fluid, I is the analogue of the electron-electron scattering length in traditional 
solid-state physics. Without quasiparticles, it is best interpreted as the minimal length scale at which a 
hydrodynamic description is sensible. The requirement that p vary slowly is often written as 


dxP 


1 

« 7 ’ 


( 7 ) 


though this should not be taken literally {p may vary slowly through p = 0). The requirement we will 
assume henceforth in calculations is that, in Fourier space, /u(k) is only non-negligible for |k|^ < 1. It is 
not necessary that p be approximately the same at all points at space:^ disorder can be non-perturbative, 
with hydrodynamic coefficients such as viscosity and charge density, contained within and in (6), 
varying substantially over distances large compared to see Figure 2. This was noted in [11] as well. 

In a quantum critical theory of dynamical exponent z, one finds 

I ~ ( 8 ) 

®|/i(xi) — /i(x2)| can be comparable to, or larger than, |p(xi)|, so long as |xi — X2I S> I- 
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by dimensional analysis [69]. Note that (7) does not hold as T —t- 0 - it is thus important that we 
are considering the finite temperature response of the QFT. (8) may be modified in models where the 
hydrodynamic limit can persist in regimes where p, ^ T, such as in holography [70] (in this particular 
model one seems to find I ~ or when the expectation value of neutral scalar fields is large. Explicit 

computations of I are beyond the scope of this paper but are necessary to properly understand the regime 
of validity of hydrodynamics. A conservative requirement is certainly to fix the background temperature 
T to be uniform and large enough that ^ S> T~^/^ at all points in space, but this may be too strict, as 
we will see in holographic models. 

More liberally, one could only require that ^ ^ I hold locally, with short wavelength disorder in “hot” 
regions of space with small I, and long wavelength disorder in “cold” regions of space with large 1. So 
long as the solution to the hydrodynamic equations of motion itself varies slowly in the cold regions with 
large I, then our hydrodynamic formalism should be an acceptable description of transport. 

When (7) holds, it is a sensible approximation (and standard in condensed matter physics) to assume 
that thermodynamic and hydrodynamic coefficients, such as viscosity rj or charge density Q, are local and 
depend only on /x(x). We then - very crudely speaking - put together pieces of homogeneous fluid of 
width whose equations of state are translation invariant, and smooth over the fluctuations from piece 
to piece. Our approach to transport is to focus on the response of the low energy, hydrodynamic degrees 
of freedom, exactly treating their evolution across the slowly varying background fluid, as we now detail. 

For holographic theories, we do not need to make the assumption ^ ^ I, or the assumption that all 
transport coefficients are functions of p alone. It is remarkable that the mathematical framework we 
develop in this paper is nonetheless applicable to so many holographic computations. 

2.1 I Linear Response: A Warm-Up 

Let us begin with some simple calculations to get an intuitive feeling for hydrodynamic transport. We 
work with first order hydrodynamics, and will justify this later in the section. We will also assume that 
our theory is isotropic, another assumption which we relax later. 

A first simple case to consider is when the only slow dynamics in the system are of charge. As we 
mentioned previously, in this case the dynamics reduce to the solution of a diffusion equation: 

dibJi = di ((t1°"(6E, - 5,5/2)) = 0, (9) 

where bEi is the infinitesimal, constant, externally applied electric field, 6 J, is the infinitesimal electric 
current, and bp is the infinitesimal local chemical potential, excited in response to bEj. (7^°^(x) is a 
transport coefficient which is inhomogeneous in space, and can be interpreted as the local conductivity of 
the theory. This approximation is well known in condensed matter physics, and we mentioned it in the 
introduction. Note that chemical potential gradients are equivalent to electric fields in the hydrodynamic 
equations of motion. The electrical conductivity of such a system is defined as 

E[bJi] = aijbEj, (10) 

where 5 Jj is evaluated on the unique solution to (9) where 6/i obeys sensible boundary conditions (e.g., 
periodicity when disorder is periodic with period L), and we have denoted with E[- • •] a uniform spatial 
average. If the disorder is isotropic, then atj = abij. Note that a is not equivalent to E[(t^°^]. Henceforth, 
we will define 

bfi = bp — XjbEj, (11) 

so that (9) can be written compactly as 

dibE = -di (a'°=5,5/r) . (12) 
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Let us now account for convective transport - this means that momentum is a long lived quantity 
and must be included in hydrodynamics. If we neglect thermal transport, then we must modify (9) to 
account for the convective contributions to charge: 


dibJi = di{Qbvi - a^dibn) = 0, (13) 

where Q is the charge density, and bvi is the velocity field in the fluid, ^ is a “quantum critical” 
transport coefficient corresponding to the flow of a current in the absence of any velocity field [ 8 ]. The 
momentum conservation equation allows us to determine and we will show more carefully below that 
this equation is the following analogue of the Navier-Stokes equation: 


Qdib^i = dj 


(^djbvi + rjdibvj + bij 



(14) 


with 77 the shear viscosity and C the bulk viscosity. As before, Q, ( and r] can all depend on position 
X, though not in an arbitrary way. As we discussed previously, the x-dependence of all these coefficients 
is hxed by their dependence on as determined in a locally homogeneous fluid. 


2.2 I Linear Response: Complete Theory 

Let us now describe the complete linear response theory which includes the response of temperature, 
chemical potential and velocity to external fields. 

About our background fluid we perturb the system with an infinitesimal electric field 


bEi = E [-^^b^J] , 


and an infinitesimal temperature gradient 


5C* = E 



Defining the heat current 




we find that the heat current is conserved (divergenceless): 


5*g* = diT^^ - diifir) = -FuJr - JAn = o. 


The thermoelectric response of the theory is given by the matrices: 

( E[6J,] \ ^ / aij aijT \ ( bE^ \ 

V E[5g*] ) V ) V 50 ) ■ 

Let us define 


5 ^" 


\ 

T-^bT ) ’ 
bEi \ 

J’ 

bJi \ 
bQi J ’ 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 a) 

( 20 b) 

( 20 c) 
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where the a vector index denotes charge (q) or heat (h). Note that we may write bold-face vectors below, 
but this always refers to spatial indices only - we will always write out the a index explicitly in equations. 
We then may write (19) as 

E[6 (21) 

We write down the gradient expansion of hydrodynamics to first order in derivatives acting on 6T, 6/i 
and 5Vi, by expanding stress tensor and charge current in terms of the linear response 6T, 6/i and 5vi of 
the fluid. The charge and heat conservation equations of the fluid may be written as 


0 = di (^Q5vi - crfjdjbn - afjdjbT^ , (22a) 

0 = di (rSbvi - Tafjdjbfi - , (22b) 

which we henceforth package into the more compact form 


0 = j)" 


di 





(23) 


where 



(24) 


where Q is the electric charge density and S is the entropy density. 27“-^ correspond to diffusive transport 
coefficients that couple charge and heat flows to gradients in 6^ and 6T, even in the absence of any 
convective (non-vanishing Vi) fluid motion. In particular, = afj corresponds to “quantum critical” 
conductivity and is typically assumed to vanish in a non-relativistic theory without particle-hole symmetry, 
as in [11]. = afj corresponds to an intrinsic diffusive conductivity that couples charge and heat flows. 

In standard non-relativistic theories, only Efj^ = Rfj is nonvanishing, as in [11]. All three are non-zero in 
relativistic systems [8] . We assume that 




(25) 


and that locally 27 be a positive definite matrix (note ai and /3j group together for purposes of matrix 
inversion). This is sensible from the point of view of the second law of thermodynamics. Indeed, in 
isotropic theories, the second law provides more constraints on these transport coefficients than (25) 
alone [8], but we can and will relax these constraints in the technical formalism we develop without 
substantially altering any physical content. This loose treatment of entropic constraints proves useful in 
Section 4. 

The momentum conservation equation becomes® 


dibP + djbTij = b [SdiT Qdifi] + djbTij = dib (^"p") - dj [rjijkidibvk] =bQdifi, (26) 

where Tij is the viscous stress tensor, P is the pressure, rjijki is the viscosity tensor with symmetries 

Vijkl — Vjikl — Vijlk — Vklij: 

and we have used the fact that thermodynamic relations imply 

diP = Sd^T + Qdi^i, (28) 


Now since T is constant on the background, and as the background p is simply given by p, we cancel the 
two 6Q terms,and we are left with 


0 = - dj [rjijkidkbvi] . (29) 

®Note that &Q = {dQ/dv)bii + {dQ/dT)hT. 
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In the above equations, p, S and tj are all smooth functions of p(x), varying on large length scales 
compared to 1. 

(28), along with (6) and the fact that = Tij = ^ on the background solution, demonstrates that the 
background solution to the hydrodynamic equations indeed exists. = 0 on the background because 
the hydrodynamic equations only depend on p — /i, which identically vanishes [8]. Though we expect 
that disorder implies that p", 27“-^ and r}ijki are all functions of p alone, we will not comment further the 
precise nature of this dependence, and a microscopic computation is necessary in general. 

(23) and (29) are linear and have a unique solution when subject to appropriate boundary conditions. 
These boundary conditions will be periodic boundary conditions in a large box of size L in all directions, 
up to non-trivial gradients 6T" = E[—We also stress that only enters the equations of motion 
through derivatives - this is crucial in order for the linear response problem to be well posed on spaces 
that are periodic or compact. Henceforth we will drop the 6 so as to avoid clutter, with a few exceptions. 

Having imposed these boundary conditions, we prove in Appendix A that, for any hydrodynamic 
transport computation. 


olB (3a 
% = ■ 


(30) 


This is referred to as Onsager reciprocity, and is a non-trivial consistency check on this framework. Note 
that this condition is violated when time-reversal symmetry (in the microscopic Hamiltonian H) is broken, 
e.g. by a background magnetic field [8]. We do not consider this possibility in this paper. 

As mentioned previously, we have truncated the hydrodynamic gradient expansion at first order. Let 
us give some sensible, though non-rigorous, justifications for this. The hydrodynamic gradient expansion 
can be organized as follows: 


% 






Jf - p“vi ~ + 


r(2)a 


(31a) 

(31b) 


corresponds to the coefficient of the stress tensor carrying n spatial derivatives; similarly for 
This is a qualitative statement - the basic idea is that the coefficients of Tlj ~ Pe/u at order in 
derivatives, e.g., with e the energy density and v a velocity scale such as the speed of sound, and so we 
have extracted out the overall scaling in I above. Assuming that the solution ^ and Vi varies over the 
length scale ^ ^ I, we see that higher derivative corrections to the charge, heat and momentum currents 
are suppressed by powers of l/^, and thus can be neglected. In the special case where diffusive charge 
and heat transport dominates, this argument can be made rigorous. When the convective contributions 
cannot be ignored, this argument is not rigorous - not all coefficients /?, E and rj scale as the same power 
of I, in general, and so it is not obvious to prove that and vt must vary on the length scale However, 
this is still a plausible assumption - any rapidly oscillatory ^ and Vi on short length scales compared 
to ^ seems non-sensible as a static solution, since static solutions to dissipative hydrodynamics tend to 
be “as close as possible” to equilibrium, given the boundary conditions; the variational methods we will 
develop in this paper also suggest that it is unlikely to have fast variations of ^ and Uj. This general 
framework readily generalizes to account for higher derivative corrections to hydrodynamics, if one wishes 
to directly include them, but we will not include them in this paper. (23) and (29) are not well-posed 
until we include first order corrections to hydrodynamics, so it is necessary to work at least to this order 
in the gradient expansion. 

In the absence of other dynamical sectors of the theory, it is necessary that either 5 or Q be position 
dependent in order to obtain finite thermoelectric conductivities. Indeed, if both S and Q are constants, 
there is a zero mode in (23) and (29) corresponding to uniform shifts in Vi and 77)“. This zero mode is 
responsible for infinite dc transport coefficients in a fully translation invariant theory. Mathematically, we 
could break translation invariance only in A or t/ and still have this zero mode. However, in a microscopic 
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theory rj, S and Q are not arbitrary but are fixed by equations of state that relate these parameters to 
ft, so in general both will be inhomogeneous. Alternatively, as we will discuss in Section 3.1, it is possible 
to add other dynamical disordered sectors of the theory which lead to finite conductivities even when S 
and Q are constants. 

Let us also briefly mention the issue of momentum relaxation times. In many holographic mean field 
models of disorder, the momentum relaxation time can be parametrically fast [53]. In these hydrodynamic 
models, the “momentum relaxation time” is parametrically slower than the mean free time (l/T in most 
quantum critical models). We do not explicitly compute this momentum relaxation time, and a single 
momentum relaxation time will not be easily definable when disorder is non-perturbative. We will see 
that it is possible to spoil coherent transport despite parametrically long lived local momentum currents. 

It is also worth stressing that henceforth, when we refer to “hydrodynamic” transport we refer to the 
transport equations being written in terms of (23) and (29). Remarkably, essentially all of our results 
rely only on the structure of these equations being obeyed, and not on S and Q obeying thermodynamic 
Maxwell relations - holographic horizon fluids do not obey any obvious Maxwell relations. If one has a 
microscopic system of interest, with known equations of state, they may simply take the above equations 
and numerically solve them. So the point of this paper is less to describe complicated (numerical) 
solutions to these equations of motion, but rather to elucidate simple and universal physical consequences 
of hydrodynamic transport: first through exact results for weakly disordered theories, and then through 
a combination of rigorous bounds and heuristic arguments for strongly disordered theories. Numerical 
solutions to these equations, and a discussion of their relevance to realistic quantum critical systems, will 
be presented elsewhere. 



Weak Disorder Limit 


In this section, we specialize to the weak disorder limit in which slow momentum relaxation dominates the 
conductivities. In this limit we can make direct contact with the memory matrix formalism [13, 14, 15], 
and provide physically transparent derivations of many previous results derived within this formalism, 
in the overlapping regime of validity of hydrodynamics and the memory function formalism. We simply 
quote the results of this approach (see e.g. [56]): if the Hamiltonian of our weakly disordered system is 

H = Hq- / d'^x /i(x)0(x). (32) 


with Hq translation invariant, and O an operator in the theory coupled to the field h, then the memory 
matrix formalism predicts, at leading order in perturbation theory: 

r _ _-i-i 


ol8 


E [p“] E 





lm{Goo{'k,uj)) 

(jj^O u 


= E [p“] E 




(33) 


L k 


In some models of strange metals appropriate for real world modeling, some care is required in defining Q 
[18]. Formally, the memory matrix formalism is exact, but the formalism does not appear to be tractable 
in practice beyond leading order, in higher dimensional models. 

Let us briefly note our conventions in Fourier space. Fourier transforms are defined as 


0(k) = L y'd'‘xe-‘'‘-=‘0(x). 


(34) 


We will often assume that disordered sectors of the fluid have (zero-mean) Gaussian fluctuations: e.g., 
quenched disorder on O would scale as 


Ed[0(k)] = 0, 


(35a) 
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(35b) 


Ed[0(k)0(q)] = ^5k,_q, (|k|^<l) 

where ^ 1 represents the number of Fourier modes which “meaningfully contribute” to 0(x), and Ed 
denotes an average over quenched disorder: 


N ~ 



(36) 


These definitions are chosen so that E[C)(x)^] = Vq. We will use these conventions throughout the paper. 
Such disorder is consistent with (7). At a typical point in the fluid, 


d^O 2 E[{d,Of] 1 

O ^ E[C>2] "" yj ^ 

To obtain the numerator in the third step above, it is helpful to go to Fourier space, and note that 
|k| < for all non-negligible modes. 

In [8] and many subsequent works, within the hydrodynamic approach to transport, the momentum 
transport equation is modified to 

= — + ■■■, (38) 

where T** is the momentum density, and r is a phenomenological relaxation time that is subsequently 
computed using memory functions.' However, as we will see in this section, at least for dc transport, it 
is actually not necessary to add in r by hand. With weak disorder, the dc transport can be accounted for 
exactly from hydrodynamics, and (33) recovered provided that the equations of motion properly account 
for disorder. 

Interestingly, our hydrodynamic approach requires the disorder is always long wavelength compared 
to I, and the memory matrix formalism does not require this restriction. Nonetheless, at leading order 
in perturbation theory, we will recover the exact memory matrix formula for transport coefficients from 
hydrodynamic considerations. It is also worth noting that the memory function approach is equivalent 
to holographic computations of transport in their overlapping regime of validity [58]. Thus, all three 
approaches give the same picture of transport, which is best physically understood in terms of this simple 
hydrodynamic framework (notwithstanding the regime of validity). 



3.1 I Disorder Sourced by Scalar Operators 

Let us begin with the case where the operator O is a scalar field. We assume that all hydrodynamic 
coefficients are x-independent - h is the only disordered parameter. In this case, (29) must be modified: 

+ diTij = 6 Odj h. (39) 

We place a 5 on O to distinguish the response due to the electric field from the background. The scalar’s 
static equation of motion is [7] 


y'd'^y G(p^(x-y,a; = 0)C>(y) =/i(x). (40) 

^In contrast, [56] has recently used the memory matrix formalism to recover hydrodynamic transport, with expressions 
for the phenomenological r and other thermodynamic coefficients expressed in terms of microscopic Green’s functions. This 
is a more complete approach to the problem, but does not generalize easily to higher orders in perturbation theory. 
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The Green’s function Gqo is the retarded Green’s function of the translationally invariant Hamiltonian 
Hq-. in position space, 

Goo(x, t) = i0(t)([O(x, t), 0(0,0)]) (41) 

with the average (• • •) taken over quantum and thermal fluctuations, and 0 the Heaviside step function. 
We emphasize that while Gqq is the true quantum Green’s function of O, and an intricate quantum 
mechanical computation may be necessary to compute it, Gqq does play the role of the coefficient of 
proportionality in the linear response of macroscopic, thermal expectation values. 

Let us make an ansatz for the solution to the hydrodynamic equations, and show that it is consistent 
with all conservation laws. Gur ansatz is that the only divergent terms in linear response, as h —)• 0, 
are v ~ h~‘^ and 60 ~ h~^. Furthermore, at leading order v = vq is a constant. All other spatially 
dependent response is 0{h?) and we will see that it can be neglected in the computation of aij at leading 
order. 

The leading order response of O in the h —?■ 0 limit is best computed in the rest frame of the fluid, 
which has shifted. It is simplest to Fourier transform to momentum space as well: 


0(k, k • Vo)co—moving — GQQ(li, k • Vo)/l(k, k • Vo)co—moving- 


(42) 


Here, everything is measured in the co-moving frame of the fluid, and so the only non-vanishing h (and 
therefore O) will have this special relation between k and w. Note, of course, that 0(k, — k-vo)co-moving = 
0(k), as measured in the original rest frame, and similarly for h(k). We wish to keep only the linear 
response coefficient, proportional to vq:® 


50(k) 


-h{k) 


dGooiKO), 

du 


■ Vo 


—i/i(k) 


Im(Goo(k,a;)) 

uj—>-0 UJ 


k ■ Vq. 


(43) 


In the latter equality we have used reality propeties of Green’s functions, and assumed analyticity near 
the real axis for k 7 ^ 0 . 

Now, let us study the momentum conservation equation, averaged over space, so the derivatives of the 
stress tensor do not contribute: 


0 = ^50(k)(-%h(-k)) 

k 

Im (Goo(k,a;)) 


^ kikjlh(k)f^ 


lim 


CJ 


voi + = -Fijvoj + 


(44) 


At leading order, the electric current is uniform: 

jr ~ p"vo^ = P^P^rr^Ff, (45) 

which gives us (33). It is straightforward to generalize to the case where there are multiple types of scalar 
fields. 

Let us now argue that the ansatz (and thus results) we have found are self-consistent. If we do not 
average the momentum conservation equation over space, then the bOdjh term is not translationally 
invariant, and this will induce corrections to T, p and v which are spatially varying. However, 50 ~ h~^ 
and so these spatially varying corrections will be ~ /i*^. Indeed, it is easy to see that (23) and (39) are 
consistent with the leading order inhomogeneous response (except in 60) arising at this subleading order. 
Thus, our ansatz is indeed correct in the asymptotic limit h —)• 0, and we have derived from hydrodynamic 
principles the momentum relaxation times derived via the memory function formalism. The computation 
above is completely analogous to the holographic computation of [58]. 

®Note that this linear term in perturbations is parametrically large in h- but it is linear in • Our perturbative parameter 
is first and foremost since we are computing a linear response transport coefficient. And while the background may also 
be treated perturbatively in disorder strength, we must take —>■ 0 before /i —>■ 0. 
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3.2 I Disorder Sourced by Chemical Potential 

In this section, we consider the case where 


Ai = Aio + eA: 


( 46 ) 


with E[A] = 0, and e <C 1 a small perturbative parameter. Alternatively, we can write 


(47) 

with = T5 and = 0 - this will be more compact notation for subsequent manipulations. We will 
denote /5“ with the fluctuations in the charge and entropy densities associated with A ~ hr general both 
will be non-zero. For simplicity, we assume that the background fluid is isotropic, though the technique 
certainly generalizes (but with more cumbersome calculations). [42] studied similar problems employing 
hydrodynamic Green’s functions in the memory matrix formalism. 

Again, we make the ansatz that the only response at 0(e“^) is a constant velocity field vq, so that 
is again approximated by (45). At 0(e“^), there are x-dependent corrections to T, /r, and v. A similar 
calculation to before gives 






A^(k) = ^ kikjp°^{-\^){m ^)"^A'^(k), 


(48) 


with 

r?' = r/ (^2 - 0 + C, (49) 

with rj the shear viscosity and C, the bulk viscosity. We provide more details in Appendix B. 

Let us briefly discuss some simplisitic limiting cases of (48) and give some analytic insight into the 
solutions - in general, the solutions will be more complicated than what we write here. 

First, let us begin with the case with rj' ^ oo and ~ 0. Suppose further p^ are Gaussian disordered 
random variables: 


Ed[A“(k)]=0, 

Ed[A"(k)A'’(q)] = -6k,-q, 


(50a) 

(50b) 


with Ed[- • • ] denoting averages over the distribution of quenched disorder modes. Then we find 


Ed [Fjj] 


^6 

Nd ^ ^ A'qq ^ dAqq 

k 


bijE 




(51) 


Fluctuations of this quantity are suppressed in the limit —)■ oo [33], so Ed[T)j] Fij. 

An alternative simple case is thermal transport with Q = 0, and 5 ~ 5o with small variations. In 
this case, we may approximately neglect the F contributions to m as ^ oo, and we find by a similar 
calculation to above: 


Ed[FiA Ri bijE 


'rj' 

J ss J’ 


(52) 


which is similar to results discussed in [42]. 

To make contact between (48) and the memory function framework, of course, we need to compute the 
retarded Green’s functions for charge and heat. Unfortunately, the Green’s functions coupling charge and 
heat flow in a relativistic hydrodynamic system are quite messy [71], and so we will prove the equivalence 
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with (33) in a more abstract manner. We proceed analogously to the previous case - as above, we have 
showed that the leading order response of the fluid that contributes to at 0(e“^) is from a constant 
shift to the velocity. 

The retarded Green’s function is defined as: 

p“(k) + 6p“(k) = G“^(k, -k • vo)^^(k) - k • VO ’ 4 ^(k). (53) 

Of course, just as in Section 3.1, we must use the boosted Green’s function to obtain the linear response 
the 0(u*^) term is the response of the background fluid, and the 0{v) term is the linear response 
contribution of interest for the computation of transport - we focus on the latter (Sp") henceforth. We 
can also relate 6p“ to by the static susceptibilities, since we are perturbing about a translationally 
invariant state: 

6p“ = (54) 

Again, we can relate vq to F" by spatially averaging (29): 


jpa 


p"(-k)ifci6^“(k) = Y / 0 “(-k)i/ci (x 

k k 


kik, 

k 


■ Im(G“/^(k,a;))' 

lim- 


SI^{k)voj. 


(55) 


It is straightforward to read off Fij from this equation, and we see that it agrees with the generalization 
of (33) to multiple disordered quantities - though of course at this point we use the fact that only / 0. 
Since Jfi ~ p°‘voij we reproduce the results of the memory matrix formalism. 

We have worked through two specific examples of deriving (33) from hydrodynamics. Of course 
one may need to generalize further, but it should be quite evident from the derivations above that the 
agreement between the memory matrix formalism and our hydrodynamic framework will persist. 

It is possible to compute the transport coefficients at higher orders in perturbation theory, where 
the memory matrix formalism has become unwieldy enough that such a calculation has not yet been 
attempted. Even at next order in perturbation theory, the corrections to the conductivity become quite 
messy. We discuss the general structure of higher order computations in Appendix B. The key point is 
that organizing the perturbative expansion in a hydrodynamic framework is straightforward in principle, 
albeit messy to carry out in practice. This framework provides quantitative predictions for memory matrix 
calculations at higher orders in e, for classes of models perturbed by fl. 


3.2.1. Breakdown of the Resistor Lattice Approximation 

Finally, let us compare the results of this subsection with the “resistor lattice” approximation: 


Ji ~ -fTij(x)5j/r(x), (56) 

with aij taken to be a local function, determined in terms of local properties of the fluid, and 

and —djpL the local electric field in the sample. Of course such a function cr may be found by solving a 
linear algebra problem at each point in space, and so the question is whether this is a useful statement 
- namely, whether crij(x) can be computed by appealing to local properties of the disordered QFT (on 
length scales large compared to 1). Essentially, can we integrate out Vi and T, and be left with a local, 
dissipative description of electrical transport in terms of p alone? 

This is impossible in the weak disorder limit, though our comments appear to have broader validity 
whenever viscosity cannot be neglected, djp is actually inhomogeneous at leading order e~^, and must 
be expressed in terms of a non-local integral over Q(x) = /9‘^(x) ~ x'^‘^e/l(x), with x*^*^ the charge-charge 
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susceptibility, assumed to be constant at this order in perturbation theory. This is derived in (101b), 
with the X and Y corrections in (101b) vanishing at leading order in e; in this equation, the leading order 
behavior of is “local in Fourier space”, and becomes non-local in position space, in terms of the original 
disorder fi. It is therefore generally impossible to find a function atj ~ e~^ expressable in terms of jl or 
its (low order) derivatives, such that Ji = aij{—dj^) = constant ~e-2 (at leading order). 


Comparing to Holography 

Let us now compare with holographic results. Many holographic results, valid in the weak disorder limit, 
are equivalent to the memory matrix results [58] - and therefore our hydrodynamic framework. So our 
focus here will be on non-perturbative holographic results. 

Our discussion of holography is brief - for further details, consult the excellent reviews [20, 21, 22]. 
Holography refers to a conjectured duality between a classical gravity theory in fi-|-2 spacetime dimensions, 
in an emergent anti-de Sitter (AdS) space, and a strongly coupled QFT in d spatial dimensions. The 
strongly coupled QFTs (in every case where we know them explicitly) are large-matrix models, and can 
be thought to “live” at the boundary of AdS. Making the gravity theory classical is equivalent to sending 
the bulk Newton’s gravitational constant to 0, and this makes bulk quantum gravity fluctuations negligible. 
This corresponds to the limit iV —)• oo in the dual theory. However, unlike vector models, these matrix 
models do not behave at all like free theories, and encode rich quantum critical dynamics. The nonlinear 
dynamics of gravity is dual to the stress-energy sector of the boundary theory. Furthermore, studying finite 
temperature dynamics becomes simply related to studying the dynamics in a black hole background. These 
black holes will be assumed to have the same planar (or toroidal) topology of the boundary theory. Adding 
a finite charge density in the boundary theory is dual to adding a bulk U(l) gauge field, and charging 
the black hole under the associated U(l) charge. Of interest for us in this paper is that holographic 
models can further be used to add strong disorder in addition to finite temperature and density. We are 
interested in modeling disorder explicitly, and so the bulk geometry becomes inhomogeneous and rugged. 
At small temperatures, the black hole gets pushed “farther back” into the emergent bulk direction, and 
the bulk fields become significantly renormalized, with higher momentum modes usually decaying away, 
as depicted in Figure 3. 

The precise duality allows us to compute correlation functions in our unknown QFT by solving grav¬ 
itational equations instead. The basic idea is as follows: correlation functions of the stress-energy tensor 
and the U(l) current are respectively related to bulk Green’s functions of the metric gMN and a 
gauge field Am, all computed at (classical) tree level. We are using MN etc. to denote all coordinates, 
including the bulk radial coordinate. For example, to compute the electrical conductivity, we add an 
explicit infinitesimal source for the bulk field A, at the AdS boundary, and then compute the expectation 
value of the current® in the field theory in the background perturbed by Aj. 

These computations are often intractable analytically. However, in the simple case of dc transport, 
many analytic computations of dc transport in holography are performed using the membrane paradigm 
[72]. In this case one can show that there is an analogous “electric current” flowing on the black hole 
horizon, which is also conserved and whose average value equals that of in the boundary (an analogous 
more complicated story holds for the heat current). The resulting computation of the conductivities then 
depend only on black hole horizon data. It is natural to conjecture that we should solve the hydrodynamic 
transport problem using the equations of state of an emergent “fluid” whose equation of state is related 
to local properties of the black hole horizon. Indeed - up to some subtleties we will see shortly - this is 
the case. 

It is remarkable that first order hydrodynamics captures the transport problem on the emergent black 
®This expectation value is also encoded near the boundary of AdS. 
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Figure 3: A qualitative sketch of holography. A finite temperature T and density boundary 
theory is dual to an emergent gravitational theory in one extra spatial dimension (depicted 
above). Strong disorder in the boundary theory (depicted in green) backreacts and leads to the 
formation of a lumpy charged black hole of Hawking temperature T. The emergent black hole 
horizon is curved and is denoted in black. The membrane paradigm suggests that dc transport 
can be computed in an emergent fluid living on the horizon, which can undergo renormalization 
relative to the “bare fluid” in the boundary theory. 


hole horizons in holography, in all nonperturbative computations to date. As we will see, however, this 
emergent horizon fluid is not locally equivalent to a fluid in the boundary QFT - instead it has undergone 
non-local renormalization, through the radial evolution of the bulk fields and geometry. In addition, 
while we previously had to make assumptions that the disorder correlation length was large to justify our 
hydrodynamic formalism, no such justification is necessary (at least, a priori) for the holographic results 
to be valid. When ^ oo in a holographic model, the differences between the horizon and boundary 
fluids become negligible, as was found in [73]. 

Let us begin with the striped models studied in [28],^^ which consider gravitational solutions to the 
AdS-Einstein-Maxwell system: 

S = Jd‘^Xy/^^R + 6 -r) 

(we have set Newton’s constant, the bulk charge, and the size of AdS equal to 1 to follow their notation). 
Above R is the Ricci scalar, and F is the Maxwell tensor associated with the bulk gauge field. In these 
models, translation invariance is only broken in the x direction, but the boundary theory lives in d = 2. 
Let us summarize their results briefly. They found that (after inverting their thermoelectric conductivity 
matrix): 


TKxx 
T O-xx 

^XX 


= aS = [b®] , 

= at = 


A 


= 


^ Hi 




tt_ 

Hu S 




(58a) 

(58b) 

(58c) 


'^°Note that their results simplify, in some special cases, to analytic results derived in [26, 34]. The case where Q = 0 was 
also studied in [30]. 
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where 


X = E 


B 





E [e^] - E 


R_^ 

Hu 


(59) 


and B, at, Hu and S are data associated with the solution of classical Einstein-Maxwell gravity, near the 
horizon of a black hole, as detailed below. The near horizon geometry in their coordinate system was 


ds^ 


Htt{x) 

HiTr 


dr^ — A'KTrHu{x)dt^ + S{x) 




+ e' 


-B(x)dy2 


(60) 


with r the radial coordinate, and r = 0 denoting the location of the black hole. The bulk gauge field only 
has a time-like component, whose value is at- T denotes the temperature of the boundary theory. 

Let us postulate the following equations of state for the emergent fluid on the horizon: 


s 


(61a) 


^ Hu 

(61b) 

= 1 , 

(61c) 

^ ^hh ^ Q 

(61d) 


The first of these equations is the canonical result for r]/S in a strongly coupled theory [74], which also 
holds for the charged black holes here [75, 76], in the translationally invariant limit, though this universal 
ratio can be different in mean-field disordered black holes [53]. The last of these equations was argued 
to occur in holographic models in [56], by matching ui = 0 results of massive gravity. More recently, it 
has been pointed out that this is not the correct interpretation of in the boundary theory, and this 
becomes discernable at finite u [54]. However, we will see that this prescription can describe correctly an 
emergent fluid, associated with data on the horizon, whose hydrodynamic response is equivalent to (58). 
The inequivalence of the boundary fluid and this emergent “horizon fluid” is an important subtlety, and 
one we will not resolve in this paper. 

We also need to make two more assumptions. The first is rather simple - let us suppose (61) is valid 
for the disordered model, with x-dependence trivially put in; e.g., r]{x) = S{x). This is in accordance 
with our logic in Section 2. The second assumption is that the boundary fluid lives on a curved space 
with metric 

ds^ = 7 ijdxMx-^ = e^^^^dx^ -|- (62) 

which differs from (60) by a conformal rescaling. Intuitively this can be argued for on the grounds that 
S determines S and rj, and therefore should not determine the boundary metric jij, since we expect that 
'jij = bij in a translationally invariant isotropic model, even though 5 7 ^ 1 in general. We compute the 
thermoelectric conductivties for such a fluid using special techniques for striped systems, discussed in 
Appendix D, and we hnd 


' XX 

-i\qh 




[cr 


J XX 


1 


V J XX 


dirT 

1 

(47rr)' 


-E 


O't 

^ H 
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rE 
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Ht, 


~rr~ 1 + v; 1 dxB — 


S 


dx^S 

S 


Inverting this matrix returns (58), the exact result found holographically in [28]. 


(63a) 

(63b) 

(63c) 
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More recently [29] has generalized the results of [28] to include the effects of a dynamical scalar field 
in the dual theory. In general this scalar must be consistently included within hydrodynamics, and so we 
must consider a more general theory to connect with these results in the generic case. 

In the case where translational symmetry is broken in multiple directions, there is an important 
subtlety. It turns out that the local “current” in the emergent horizon fluid is not equivalent to (J) in 
the boundary theory. This is true only after taking a spatial average [60]. To relate the “current” in the 
horizon fluid to the current the boundary fluid, one must add a non-local integral over the bulk direction. 

It was recently shown [60] in a more general context that dc transport in holography reduces to solving 
“hydrodynamic” equations on the black hole horizon. [60] interpreted the resulting fluid equations as an 
incompressible Navier-Stokes equation. [35] points out that these equations can also be interpreted in the 
framework of the present paper. 

These examples suggest that - while the hydrodynamic framework of this paper is extremely helpful 
providing some physical intuition to these non-perturbative holographic results - this story is not complete. 
Importantly, however, much of the variational technology that we develop can be directly applied to 
holographic models. 



Strong Disorder 


We cannot be as rigorous in the strong disorder limit and give closed form expressions for the conductivity 
matrix. Nonetheless, we will develop simple but powerful variational methods that allow us to get a flavor 
of transport at strong disorder, by providing lower and upper bounds on the conductivity matrix. 

We focus on the discussion of in an isotropic theory in this section. However, the techniques 
developed below may be used to compute all thermoelectric transport coefficients. Our discussion is 
therefore not an exhaustion of all possible physics contained in the hydrodynamic formalism, but simply 
a demonstration of what we believe is a general feature of hydrodynamic transport: a crossover from 
coherent (Drude) physics to incoherent behavior as disorder strength increases. 

We present the mathematical formalism in the subsections below - explicit examples of calculations 
may be found in Appendix C. 


5.1 I Power Dissipated 

Define “voltage drops” of each conserved quantity in each direction as 

I/“ = = 0, x_Li) - = L, x_Li)- (64) 

Recall our boundary conditions on it must be periodic up to linear terms. We also define the net 
currents flowing in the i direction via 


= 


d'^-^x 


fixed Xi 


(65) 


Note that by current conservation. If can be evaluated at any x*. Since V and I are determined by the 
solution to a linear response problem, we can relate them via a conductance matrix (do not confuse this 
G'"4 with the retarded Green’s function defined previously) 


If ^ Gffvf, 


( 66 ) 


or via its inverse, the resistance matrix 


'^ 1 - bj ■ 


(67) 
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is by definition related to the dc transport coefficients: 

,^«/3 _ Td-2aP 
^ ^ij ■ 

We claim that the power dissipated in the system is simply given by 

V = 


( 68 ) 


(69) 


Let us verify this. Energy is dissipated^^ locally via the dissipative {U and tj) terms in hydrodynamics: 


V= / d'^x + mjkidjv, 


(70) 


We integrate by parts on the second term and use (29) (recall Vi obeys periodic boundary conditions): 


V = 


j d'^x = - y d'^x 


But 77)“ is a conserved current, and so 

V = - j d'^“^x <P'^niX = = 0) - ^“(xi = L)) = 

i 

with Ui the outward pointing normal. 


(71) 


(72) 


5.2 I Lower Bounds 

Let us begin by discussing the lower bounds on conductivities. These are by far the more important 
bounds to obtain, because - as we will see - they allow us to rule out insulating behavior in a wide variety 
of strongly disordered hydrodynamic systems. 

We obtain lower bounds on conductivities analogous to how one obtains upper bounds on the resistance 
of a disordered resistor network, via Thomson’s principle [77]. Similar approaches are also used in kinetic 
theory [78]. Thomson’s principle states that if we run any set of “trial” currents through a resistor network, 
subject to appropriate boundary conditions, then we can upper bound the inverse conductivity by simply 
computing the power dissipated by our trial currents. The power dissipated in the resistor network is 
minimal on the true distribution of currents, which is compatible with Ohm’s Law and a singly-valued 
voltage function. We will see that remarkably, this simple approach immediately generalizes. 

Let us propose a trial set of charge and heat currents, 77)“, which are periodic functions, and exactly 
conserved: 

9*77“ = 0. (73) 

In general, this trial function will not be compatible with a single-valued (well-dehned) ^“. We write 


rra _ /ja I /^Q; 

Ui Ui \ Ui 


(74) 


with overbars denoting the true solution of the hydrodynamic equations subject to our boundary con¬ 
ditions, tildes denoting our trial “guesses” at the true solution, and hats denoting the deivations, on all 
variables henceforth. We also impose 


d'^-^x77)“(x* = 0,L) = 0, 


(75) 


course, this would lead to temperature growth at second order in perturbation theory, so that the energy conservation 
equation (up to external sources) exactly holds at all orders. 
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as there is a true solution with the same net currents /“ as our trial We also propose a completely 
arbitrary periodic velocity field vi. Define 


V= / d'^x 


{jt - {jj - ) + rjijkidjVidiVk 


{E-^) 


(76) 


which, on the true solution, is analogous to (70). We define V (the true power dissipated) and V analo¬ 
gously. Recall V, V, "P > 0, and expand out V: 


r = 'P + 'P+ 2 / d'^x 


{E ^)")^ - p^Vi] - p^Vj] -h PijkidjVidiVk 


0 _ 


= V + V + 21C 


(77) 


Now 

(jt - P^v^) (jf - /u,) = - (E-^yf (jr - p^Vi) E^Ok^^ = - P^h) (78) 

and so we obtain, integrating by parts: 


/C 



di<P°‘J°‘ -b -b PijkidjVidiVk 

x^“ni7^“ -b J d'^x - djipijkidiVk)] Vi = 0. 


(79) 


The first term vanishes since 77)“ is periodic, and the constant gradient terms vanish due to (75); the 
second by (29). We conclude that V >P. If we define V = Ryifl^, then we obtain 

irnyy > m 


In particular, we can immediately obtain bounds for all diagonal entries of ■ Suppose that we 
have a large, isotropic disordered metal, in which case we find and = BR^hij. Then 

ZJ J IJ J 

we have the generic bounds 


1 1 

^3- < - 

Rm ~ Rqq 


1 

Rhh 


< 


G^hG-qq _ (Ghq)2 
G^h 


< G<^^ 


It is also straightforward to convert from into (a : 


(81a) 

(81b) 


V = IfR^y = [77)“ 


\RPfE 


jP 


= L'^E [77“] (a-i)“/E 


JP 

'P 3 


(82) 


To obtain bounds on off-diagonal elements is slightly more subtle. Information can be found about off- 
diagonal elements by studying linear combinations of various components of but it is not as clearcut 
as (81). 


5.3 I Upper Bounds 

Obtaining upper bounds on conductivities is in principle more simple, but quite a bit more subtle in 
practice. Let us write 

p = I d'^x + pykiTijTk) = vtcyvy (83) 
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where Tij = —rjijkidiVk is the viscous stress tensor, which has d{d + l)/2 independent components; 77 “^ 
is a matrix inverse with the first two indices grouped together and the last two indices grouped together 
(but only in symmetric combinations). It is possible that ry may not be invertible, but in this case it 
is straightforward to regulate the zero eigenvalue with an inhnitesimal positive eigenvalue and then take 
the inverse.To compute the conductance, we need to demand (as stated previously) that obeys 
^«(xi = 0) - = L) = y". 

We are going to guess a single valued trial function with the exact solution as 

before, and ^ a periodic function (recall that should be periodic up to the linear gradient terms). We 
will also guess a trial Tij = Tj + Tj, which must be a symmetric tensor, and a periodic function. We do 
not require that Tj be expressable in terms of a velocity function, or that diJ^ = 0. Let us verify the 
circumstances under which we can nevertheless hnd V > V, as before. We hnd that V = V + V + 2/C, 
with 

/C = y d'^X + V^hfijTkl) = I d'^x ((p“h, - Jr) dj" - %llTjVklmndnVm) 

= j d'^x + dj%) Vi- ^ d-^-^x mJrS^ (84) 


We have used the periodicity of Tj to integrate that term in JC by parts. The first term vanishes if we 
require that 

p»di§'^ + djfij = 0 (85) 

of all perturbations. The second term vanishes because and J^ are periodic, with UiJ^ taking opposite 
signs on each face. 

And since the integrand in (83) is positive semi-definite, "P > 0. We conclude that this forms the basis 
of a variational principle for the computation of If we dehne 


p = 


( 86 ) 


then we obtain bounds on G^^, analogously to G'^f. As before, diagonal elements of can be straight¬ 
forwardly upper bounded, and off-diagonal elements require more care. 


5.4 I Discussion of Variational Results 


Here we present a summary of the calculations performed in Appendix C. For simplicity, let E[Q] = Qq 
and Var[Q] = E[Q^] — E[Q]^ = T. The electrical conductivity will, in a disordered isotropic fluid without 
parametrically large fluctuations in 27 or 77 , be bounded from above and below by the following schematic 
bounds: 


Qn Qn 

CrQl(M) (Tq2{u)^ < O' < 0-q3(u) -h fTQ4(u)^ 


C'Qo 


I TT ^ 


V? rii{u)v? miT VsT) 


(87) 


with each cjq and 77 factor above related to “typical” behavior of and pijki respectively. In particular 
the upper bounds are quite subtle (see (126)), and so each (Tqi^ 2 , 3,4 may have complicated u dependence 
for u > Qo, and we have written (87) as we did to emphasize qualitative behavior, as discussed after (2). 

(87) proves that there is a crossover at u ~ Qo between an coherent regime when u <C Qo (translational 
symmetry is weakly broken) and an incoherent regime when u P> Qo (translational symmetry is strongly 
broken), as depicted in Figure 1. As discussed in the introduction, this is the physics found by mean held 
holographic models, and demonstration of this without a mean held treatment of disorder is a primary 
quantitative result of this paper. 
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For example, this may correspond in a conformal fluid to deforming the equation of state with a non-zero bulk viscosity. 
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Many of the statements which lead to (87) can be made quite rigorously. The lower bounds on 
conductivity are derived carefully and will be valid in a wide variety of theories. The upper bounds which 
we derive are much more challenging to evaluate analytically when viscosity is not neglected, and so we 
have made heuristic arguments to understand the qualitative physics that are non-rigorous and may break 
down in some cases. Theories with very large fluctuations in U, rj, or /?“ could render the upper and lower 
bounds to be far enough apart (perhaps parametrically so) for (87) to not be useful. Still, we propose 
that the coherent-to-incoherent crossover described by (87) is generic, and will provide some intuition 
into why this occurs. 

In (2), we ignored viscous effects, and in Figure 1, we depicted a saturating at a* when u ^ Qq. (87) 
generically confirms this picture, with cjqi < ct* < Uqs, so long as 

( 88 ) 


This inequality may or may not be satisfied, and determines whether transport may become sensitive to 
viscous effects. For example, in a strongly interacting quantum critical system of dynamical exponent z, 
we expect cjq ~ r? ~ 5 ~ [ 74 ]^ and ^ The requirement that ( 88 ) is violated is 

Tid-2}lz Xl/z^ 

When /i < T (the regime of validity^^ of the hydrodynamic approach in a typical quantum critical model), 
then it is reasonable to expect that m < 5, as most of the entropy will be associated with a background 
charge neutral plasma, and not with the deformation by a chemical potential. Rearranging (89) we find 





T{d-2)lz rpl/z^- 


(89) 


1 > - > 
r\_/ ^ t-N_/ 


CFq 


T(d-2)/z 


(90) 


Recalling that ^ 1, this sequence of equalities is satisfied for disorder on the longest wavelengths. 

However, if u <C 5 then it may be possible to have disorder on short enough wavelengths that this 
sequence of equalities is not satisfied. It is this regime where (2) and Figure 1 are valid. A viscous- 
dominated transport regime is not understood well, and a further understanding of this regime is an 
important goal for future work. 

A second assumption that went into (2) is that (A’“^)'^'^ is finite. This is true in the effective horizon 
fluid of holographic models, but need not be true in other quantum critical models. Transport in models 
where (A'“^)'*^ is infinite will be discussed elsewhere. 

Similar bounds can be found for other transport coefficients. In particular, for bounds on R, one must 
simply replace cjq with TRq, Qq with TSq = E[r5], and with T^Var[5]. It is more likely that thermal 
transport is sensitive to viscosity in a quantum critical system, as kq —)■ 0 when /x <C T [8] . 

One of the most important results we find is an exact inequality for an isotropic fluid: 


a 


> 




(no summation on a or i). 


(91) 


This can be interpreted simply as the statement that a uniform charge or heat current could flow through 
the fluid, with no convective transport, encountering this effective conductivity. And so as long as a 
current can flow everywhere locally with a finite conductivity, so can a current flow globally. (91) is 
incredibly powerful - in particular, if is strictly positive at all points in space, we have proven that 
the QFT described by this framework is a conductor. As mentioned above, the lower bounds in (87) are 

^®When /X 3> T, generally Fermi liquid theory is valid. 
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derived quite carefully, and essentially follow from generalizations of (91). Our proof that these fluids are 
conductors when is finite generalize to anisotropic theories, though the bounds become more easily 
expressed as upper bounds on the matrix 

The new approaches advocated in this paper, along with the existing mean-field literature, suggests 
that the fate of most holographic models at strong disorder - at fixed temperature T, and arbitrarily 
strong disorder - is to become an incoherent conductor, and not an interacting quantum glass. This 
is a remarkable and highly non-trivial prediction. In contrast, in metals described by (fermionic) free 
quantum field theories, there is a transition to an insulating phase at some critical disorder strength [79], 
which is zero in d < 2 [80]. Generically, interactions do lead to delocalized, conducting phases at weak 
disorder, with localization and insulating physics arising at stronger disorder strengths, in any d [81, 82]. 
It is possible that this localization transition is not observable in classical holography, which only captures 
the leading order in —)> oo, and so has taken the “coupling strength —>■ oo” limit before the “disorder 
—>■ oo” limit. 

We have always referred to these hydrodynamic models as incoherent metals. Holographic “insulators” 
discussed in the literature typically rely on scaling as a positive power of T, in a homogeneous model. 
In [29], it is likely due to stripes of such decreasing arising at low T. More generally, such insulators 
arise from the percolation of locally insulating 17'^^ regions through the effective horizon fluid.This is not 
unlike the “metal-insulator” transition of a classical disordered resistor lattice, associated with percolation 
of i? = oo resistors across the lattice [65, 83]. This is a different mechanism from Anderson localization in 
typical condensed matter systems, which is related to destructive interference of quasiparticles scattering 
off of disorder. Of course, in holographic models, the percolation phenomenon on the horizon could 
emerge from “benign” disorder on the boundary, but from the point of view of the emergent horizon fluid 
the metal-insulator transition is simply a percolation transition. We emphasize that our hydrodynamic 
formalism is still mathematically valid for dc transport in holographic insulators, due to the remarkable 
mathematical results of [60, 61]. The physical interpretation of such a fluid is an important question for 
future work, as emphasized in the previous section. There is to date no construction of a holographic 
metal-insulator transition that is unambiguously driven by (non-striped) disorder, and interpreting any 
such model in terms of hydrodynamic transport may lead to interesting insights. In simple holographic 
models, it has recently been shown that such a transition is impossible [35], and so more complicated 
models with bulk scalar fields will be necessary. 

In a non-holographic context, it is less clear whether or not our hydrodynamic formalism will be valid 
in a quantum system undergoing a metal-insulator transition, as the validity of hydrodynamics rests on the 
disorder being long wavelength. The classical “metal-insulator” transition realized by resistor networks 
[65, 83] is a crude example of this phenomenon, but relies on the only hydrodynamic degree of freedom 
being charge. 

Finally, as we are studying a strongly disordered system, it is also worthwhile to think about fluc¬ 
tuations in the transport coefficients between different realizations of the quenched disorder. As in [33], 
we expect that these fluctuations are suppressed as the size of the sample increases with L as 
with possible deviations when distributions on the random coefficients p, E and p are heavy tailed. Such 
fluctuations are classical, but this is not surprising since the dc response of our QFTs are governed by 
classical hydrodynamics. This is analogous to weakly interacting theories at finite temperature [84]. In 
contrast, a free quantum field theory has universal conductance fluctuations at T = 0 [85, 86, 87], so it 
would be interesting to ask if the T —)• 0 limit of holographic models (where hydrodynamics can still be 
a sensible approach [70]) has anomalous fluctuations in transport coefficients, in disordered models. 

the percolation problem is trivial in d = 1, the study of holographic insulators may be quite a bit richer in models 
with translational symmetry broken in multiple spatial dimensions. 
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Localization 


As we previously mentioned, many free or weakly interacting quantum systems are described by a “lo¬ 
calized” phase where transport is exponentially suppressed at low temperatures [79], at strong disorder. 
Naively, one might think that a strong coupling analogue of localization - with the associated reduction 
in transport - would exist at strong disorder. Indeed, [88] provided evidence for a possible connection in 
a holographic model. In seeming contrast to this, we have rigorously ruled out any insulating, localized 
phase in our framework (which includes many such holographic models), so long as the quantum critical 
conductivity is finite everywhere; the simple holographic models studied in the literature to date are 
described by our framework, with finite everywhere in space, in most models. 

This is consistent with known results in elastic networks and other random resistor networks. Despite 
the localization of classical eigenfunctions [89, 90, 91], diffusion and transport are possible even with 
localized eigenmodes of the linearized hydrodynamic equations. This has been shown in similar models 
without convective transport [92, 93, 94, 95]. Localization is more subtle in these systems due to the 
presence of zero modes of the hydrodynamic operators, due to exact conservation laws. Together with 
modes of arbitrarily long correlation length with finite eigenvalues, transport is possible despite classical 
localization, and so the signatures observed in [88] need not be important for dc transport. 

The finite momentum or finite frequency response of the system may be more sensitive to localiza¬ 
tion. In a simple model of disordered RC circuits, interesting new universal phenomena arise [96]. It is 
worthwhile to understand finite oj transport in the class of models in this paper as well. In particular, it 
is interesting to ask whether at strong disorder, the Drude peak found via memory matrices [56] broadens 
out enough to look “incoherent” [38], at least at small frequencies, or whether more exotic phenomena 
emerge. 

There is one other point worth making about localized eigenmodes. In a translationally invariant 
fluid, long-time tails in hydrodynamic correlation functions in d < 2 spoil hydrodynamic descriptions 
of dc transport [71]. In particular, in d = 2, the conductivity cf{uj) in an uncharged, translationally 
invariant system picks up a correction ~ log(T/|a;|), which diverges as w —)> 0 [97]. In holography, it is 
known that such long time tails are quantum bulk effects [98] , and are thus completely suppressed in the 
models described in Section 4. Since it has been argued that the memory matrix approach gives sensible 
predictions for realistic strange metal physics [17, 18, 19], and the memory matrix framework employed 
there can be interpreted hydrodynamically, one might be, a priori, concerned about whether long time 
tails can spoil dc transport in these models. If in the thermodynamic limit, all modes (except for the two 
zero modes) of the classical hydrodynamic operators are localized, as is believed in d < 2 (where long 
time tails are problematic), then the standard argument for long time tails [71] seems to fail. It would be 
interesting to explore this point further in future work. 


Conclusion 

In this paper, we have explored the consequences of hydrodynamics on the transport coefficients of a 
strongly coupled QFT, disordered on large length scales. We demonstrated that hydrodynamics can be 
used to understand the memory function computations of momentum relaxation times, which have pre¬ 
viously been derived using an abstract and opaque formalism. It is also straightforward - at least in 
principle - to compute transport coefficients at higher orders in perturbation theory, whereas memory 
function formulas only give leading order transport coefficients. Remarkably, we also demonstrated that 
many non-perturbative holographic dc transport computations can be interpreted entirely by solving a 
hydrodynamic response problem of a new emergent horizon fluid. Thus, the technology of Appendix C 
may be applied to these models, when exact solutions are not available. We still need specific microscopic 
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theories to compare with T-scaling laws in experiments, but this work provides important physical trans¬ 
parency to a large body of recent literature on transport in strongly coupled QFTs. We emphasize again 
that the memory function formalism and holography are valid in regimes where hydrodynamics should 
formally break down, and so it is strange (but useful) that hydrodynamic technology (which is readily 
understandable) can be used to help interpet these results nonetheless. 

The fact that this hydrodynamic framework can be used to interpret such a wide variety of results from 
memory function or holographic computations is suggestive of the fate of such theories at strong disorder. 
Shortly after this paper was released, it was proved in [35] that all AdS 4 -Einstein-Maxwell holographic 
models are electrical conductors, using these hydrodynamic techniques (which are valid so long as the 
bulk black hole horizon is connected). Thus, we do not expect a holographic analogue of a many-body 
localized phase [82] to exist in many strongly disordered holographic systems. Strongly disordered black 
holes have been numerically constructed recently [99, 100, 101]; hence, dc transport coefficients in these 
backgrounds, along with finite momentum or frequency response, may be numerically computable in the 
near future. 

More generally, we have also demonstrated ~ without recourse to mean field treatments of disorder, 
or to holography - a framework which generically gives rise to both a coherent metal at weak disorder, 
and an incoherent metal at strong disorder. Incoherent metallic physics has been proposed recently to 
be responsible for some of the exotic thermoelectric properties of the cuprate strange metals [102, 103]. 
One can easily imagine taking realistic scaling laws from microscopic models appropriate for cuprates 
[17, 18, 19, 57] and making quantitative scaling predictions about the strong disorder regime using insights 
from our framework. 

Hydrodynamics provides a valuable framework for interpreting more specific microscopic calculations. 
There are many natural extensions of this work: two examples are the study of hydrodynamic transport 
in disordered superfluids and superconductors, or the study of systems perturbed by further deformations 
than p,. Still, this framework has limitations. High frequency transport (in particular, ui > T) cannot 
be captured by hydrodynamics, and provides a unique opportunity for holography in particular to make 
experimentally relevant predictions about quantum critical dynamics [97, 104]. 
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Onsager Reciprocity 

In this appendix we prove that the thermoelectric conductivity matrix is symmetric. This follows 
entirely from the symmetries of the diffusive transport coefficients (25) and (27), as well as the equations 
of motion (23) and (29). 

To do this, we look for a periodic solution and Vi, up to the constant linear terms in In 
particular, we write 

= -F^Xi + ^fFf, (92a) 

Vi = (92b) 


27 


which we may always do, as the equations of motion are linear. (23) and (29) become: 


We further have 








= E 



(93a) 


(93b) 

+ 'ffSkSS + . 

(94) 


as we can now always integrate by parts inside of spatial averages. Now, let us employ (93) and (25) and 
write 




"IjP'Bk'i’l" + + dk'i’f (pTiig - s'lf ) + r; 




-,a/3 

"i-j 


(95) 


Using (25) and (27) it is straightforward to see from the previous equation that is symmetric. 



Perturbative Expansions 


Let us describe how to extend the weak disorder calculations of Section 3 to arbitrarily high orders in 
perturbation theory, in the special case where the disorder is introduced entirely through fl. This also 
gives a flavor for how to “extend the memory matrix formalism” beyond leading order in perturbation 
theory. 

Let us write ^ = hq + with e <C 1 a perturbatively small number. Within linear response, the 

fields and Vi may be written as follows: 


oo 


= 

Vi = 


E 

n=—2 


n=—l 

oo 

(96a) 


(96b) 


n=—l 


where E[v] = 0, v a constant, and single-valued. We will justify the powers of e above in our 

computation below, but for now let us emphasize that enter the computation at 

the same order. In addition, the hydrodynamic background becomes disordered: 

OO 

p” = pS + E'>“»). 

n=l 

oo 

Ef = sfbii + V , 

tj U ' / V ij (n) 

n=l 
oo 

Vijkl VjOijkl T ^ ^ 

n=l 

The background at O(e^) is translation invariant, but not at higher orders. As in the main text, we will 
assume isotropy of the leading order transport coefficeints. 

(23) and (29) may be perturbatively expanded in powers of e. We find the following equations in 
Fourier space: 

ifciP“i)(k)lli(n-i) + i^i/5of^i(n)(k) + = -A(“)(k), (98a) 


(97a) 

(97b) 

(97c) 
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i^iPo^fn)(k) + voh + kjVi(^n){^) - 2^ijkiVi(n)0^) j + Co^ijkivi(n)0^) = -^i(n)(k), (98b) 

iX;/^.^H(k)pf(-k) = -Z,(„) (98c) 


with the third equation the zero mode of the second, and 

^(-1) = 0 
ki(-i) = 0 


and, for n > 0: 


n—2 


n—1 


^{n)=^ki Y, P“n-m)(k)w*(m) +ifei ^ Pfn-m) (k " q)f 5 i(m) (q) 


m =—2 
n—1 


m=—l 




m=—l q 
n—1 

Yiin) = YY^ [i^i^fm)(q)^fn-m)(k “ q) + V{n-m){^ “ ^)kjqjVi^m){<l) 

m=—l q 


+ (^(n-m)(k - q) - ??(n-m)(k - q)j kiqjVj(^^){q) 

n—1 

Ziw = -iE E 

k m=—1 

Order by order in perturbation theory, these equations may be solved exactly: 

I = -rp Y - i (x;(k) - 




rn)(k) = -i(m(A:) ^)"^ f feiUi(„_i)p^^^(k) + X(“ )(k) - iPo^ 


. pkiY,i(jp{[Fj 


'"'”1 = (*« 


Vo - Vo hkj 


Vo 


^2 j ^■ 0 )(k)- 


(99a) 

(99b) 

(99c) 


(100a) 


(100b) 

(100c) 

(101a) 

(101b) 

(101c) 


with Fij and m“^ given by (48), with —?• p'^^y These equations have a clearly nested structure and can 

be iteratively solved. At leading order, it is readily seen that the response of the fluid is simply 


^*(—2) PoPo^ij ’ 


( 102 ) 


as claimed in the main text. We also stress that even at leading order, F°‘ and Vi are non-local functions. 

At higher orders, the A, Y and Z corrections must be systematically accounted for, and this is 
overwhelming to process by hand, especially without specific equations of state. However, this tedious 
procedure does seem easier than attempting to generalize the memory matrix formalism to higher orders 
in perturbation theory, and indeed makes predictions for such an effort. So let us at least comment on, 
qualitatively, what happens at higher orders in perturbation theory. Many terms that contribute to JP 
at higher orders in e are related to p'^^y Yf(n)^ Pin) C(n)) oi higher powers of n. As we mentioned 
in Section 3.2, we can interpret 


P“i)(k) = Re [G"‘5(k,a; = 0)] /i(k) = x“V(k). 


(103) 
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Namely, the response coefficients above are related to certain Green’s functions that can be computed in a 
microscopic model. Recall that the disorder is on such long wavelengths that we may neglect k-dependence 
in the hydrodynamic Green’s functions. So it is tempting to interpret 

/5fn)(k) = ^ ^ Re[G“‘i-‘i(ki,...,k„,w = 0)]/l(ki)---/l(k„)6k,ki+...+k„ 

~ ^ , (104) 

with G°“^ an appropriate n-point Green’s function in the microscopic theory. Similar statements may be 
made for U and r] by relating them properly to Green’s functions of and T^y, as in [71]. In the last 
step, we have used the fact that disorder is long wavelength, and so we expect Vijki{n) to 

be local functions of // in position space. 

This provides a prediction of our hydrodynamic framework which may be compared with a memory 
matrix calculation at higher orders in perturbation theory (or another method). Of course, we should 
stress that in principle, memory matrix calculations can account for corrections beyond the regime of 
validity of hydrodynamics, though in the limits we identified in Section 2, one should find that only the 
contributions described above contribute to the conductivities. 

In the above framework, it does not seem as though there are any natural cancellations between various 
terms at higher orders in perturbation theory. So this approach becomes rapidly unwieldy for computing 
transport coefficients past leading order in e. Holographic mean field phenomenology suggests that these 
corrections are all related to a single phenomenological coefficient - the Drude relaxation time r ~ in 
(1). It would be interesting to understand further under what circumstances the Green’s functions above 
undergo similar universal cancellations, and whether this is a sensible prediction of holography. 



Examples of Variational Calculations 


C.l I Upper Bounds on the Resistance Matrix 

A simple set of trial functions is 


77j“ = constant, (105a) 

Vi = 0. (105b) 


This is a guess corresponding to strong momentum relaxation, as the response of the metal is entirely in 
the diffusive sector. Employing (81) we obtain (91). 

In cases with weak disorder this bound is not strong enough to be useful, and we can do better by 
allowing for Vi to be a constant (x-independent) variational parameter. In this case, we obtain 


P{vi) = L-" + 2B^^jfvi + Cij 




(106) 


where 




= -E 
Cij = E 


(^-l)«/3 
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(107a) 

(107b) 

(107c) 












Minimizing V{vi), we find 




a 




jf<jr 


4"^ _ R“ 




(108) 


It is straightforward to see that the smallest eigenvalue of a~^ must be larger than the smallest eigen¬ 
value of a~^. A generic consequence of this result is that if the components of a~^ are not parametrically 
small in the weak disorder limit, the components of a may be parametrically large in the weak disorder 
limit. A simple example analogous to Section 3.2 is the case where A is a constant, isotropic matrix, and 
p°‘ = Pq + with Pq a constant and /5" a small perturbation with E[p"] = 0. In this case we hnd 


(a-Xf = 


(XIpo + 


(109) 


As p°‘ —>• 0, one can check that pg becomes an eigenvector of with a parametrically small eigenvalue. 
Exactly at = 0, will have an eigenvalue of oo; as discussed previously, this follows on quite general 
principles from the fact that S and Q become constant. If we invert this matrix, we find 


a 


OL^ 


PoPo 


PoPo 




+ 


( 110 ) 


To compute this eigenvalue,it is easiest to compute Pq-, and take the leading order coefficient 

in /5". The subleading corrections correspond to diffusive transport and stay hnite in the —?■ 0 limit. 

Let us compare to the exact results in the perturbative limit in Section 3.2. Technically speaking, 
we are not guaranteed that < cr“^, though this inequality is satished in this limit (assuming that 
Pq > 0). For we can write Ri pgPg/C, with C given by (48), and 


c = + /(k)<l^r(-k)(i;-'rV(k) = )). (Ill) 

k \ ' J k 


The first inequality here follows from the fact that for any vector Ui, and two positive definite matrices 
Aij and Bij, the following inequality holds: 


Ui{A + B)^Xi ^ 




( 112 ) 


To prove this, let A > 0 be a positive coefficient, and 


-^^Ui{A \B)-^Uj — —Ui{A Bki(yA -\- XB^^^Uj < 0, (fl3) 

with the latter inequality following from positive-definiteness of sums and products of positive definite 
matrices. Integrating (113) from A = 0 to 1 proves (112). 

It is also possible to find viscosity-limited bounds on the resistance matrix which can be smaller than 
the diffusion-limited bound (108), where viscosity plays no role. For simplicity, let us focus on the specihc 
case of computing thermal transport in an isotropic theory with Q = 0. (108) gives us that 


(c7"^)‘'‘"<E (A-^) 


-l\hh 


E[(^-l)hh5]2 
E[(i;-i)hh 52 ] ■ 


(114) 


compute an eigenvalue, one should first properly “re-dimensionalize” p°‘ so that all matrix elements of a°‘^ have the 
same dimension. 
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A natural guess for a viscous-dominated bound is to assume that 


Jt 


Vi 


= constant, 



TS' 


This directly leads to 

rj’ / 9*5 y 

d \TS^) 

We may employ whichever of (114) or (116) is smaller. 


(u-i)^*^ < E 


(115a) 

(115b) 


(116) 


C.2 I Upper Bounds on the Conductivity Matrix 

For simplicity, we focus on the bounding of may be bounded with an exactly analogous ansatz. 

Let us write the background charge density as 

Q = Qo + Q, (117) 

with E[Q] = Qo ^ 0 and E[Q] = 0. Let us also split into a linear term sourcing a background electric 
field, and a periodic response ip: 

= +(118) 

i 

It is easier to deal with (85) in Fourier space, so let us write (with Ei = V^jL): 

i£'*Q(k) + ^ g'*(^(q)Q(k - q) = kjTij{\s). (119) 

q 

Let us begin by assuming that ?] —)> oo, so that we may ignore the response of Eij to when computing 
(83). The only constraint we must impose in this limit is 

E[QdiS^] = 0 . ( 120 ) 

A natural guess for ip, inspired by exact results in the weak disorder limit in Section 3.2, is 

¥’(k) = (121) 

with A a positive constant constrained by (85): 

EA = E,Y.^\Q(<i)\^ ( 122 ) 

q ^ 

So far, up to neglecting viscous contributions to V, this is completely rigorous. 

For simplicity, suppose that Q(k) are disordered random variables, that are not drawn from a heavy¬ 
tailed distribution: 


Ed[Q(k)] =0, 
Ed[Q(k)Q(q)] = -6k,-q. 


(123a) 

(123b) 
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This will allow us to extract meaningful qualitative information out of bounds which will be quite opaque 
in Fourier space. A is fixed by (122) as iV —)• oo: 


A = 



Fluctuations in A are suppressed as N [33]. 
Plugging this ^ into (83) we obtain 


(124) 


XJ J 




= E 


d-2 


E 


^3 


EiEj - 2E 




Ei + E 




(125) 


Now, recall that we expect is a function of /x, and Q is a function of n as well, so let us just consider 
to be a function of Q. Then we obtain (exploiting isotropy): 


Ed[<] < ^Ed 


6 29 


d2«(e)-2^V>‘i(-k)5M+ ^ rM(-k,-k2) 


ki,k2 


di:qq(g)_2^E‘i'i(-k)® + E'^'i(-ki - ka) 

k ki,k2 


Q(ki)Q(k 2 ) (ki-k 2 

A2 

.XT. t\j<2 

Q(ki)Q(k2)’' 
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= Ed 


_ 2QoQEqq ^ dQlQ^E^^ 








(126) 


This equation holds for arbitrary functions E'i'i(Q) > 0, so long as viscosity is negligible. Our disorder- 
averaged bound on Uij is also manifestly positive, as it must be. 

In a theory with Qq —)• 0, and viscosity still negligible, our upper bound collapses to E[E'i‘^]. This 
bound may also be found at Qo = 0 by directly plugging in the satisfactory ansatz = —EiXi into (83), 
and so our bound is actually valid for all Qo- 

Now, let us consider the effects of finite viscosity. Henceforth, the discussion will be more qualitative, 
and we will not be particularly concerned with 0(1) prefactors, as it turns out to be quite difficult to 
write down a good non-perturbative analytic solution to (85). 

Let us see what happens if we simply use (121), along with a sensible ansatz for Tij. Denoting 


iF;jQ(k) ^ gi¥?(q)Q(k - q) = A(k), 

q 


(127) 


we pick in d = 1: 


in d = 2: 


and in d = 3: 


2T- _ 

f XX — ? 

kx 


^ _ 27 - _ kxAx kyAy 

I XX - - lyy - L2 , t.2 ’ 

kyAx T kxAy 

Ixy = 


1.2 I 1.2 


Txy — 


kxAx T kyAy kzAz 

‘2kX ky 
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(129a) 

(129b) 


(130a) 
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(130b) 


= 


Tyz — 




‘2kxkz 

kx*^X “1“ kyjAy “t" kz*^Z 


2kykz 


(130c) 


In all of the above cases, the equations are only valid at k = 0, and we take the zero modes to vanish. 
In d = 3, we make the stronger assumption that, e.g., Txy = 0 whenever kx = 0 or ky = 0. It is 
now straightforward to (qualitatively) see what happens. The first contribution to the conductivity is 
unchaged from (126), and the average viscous power dissipated scales as 





k^O 

V 






V 


rju^ 


ki ,k 27^0 

EiEi 


kik2 


(131) 


where we have been schematic and neglected tensor indices on rj. To obtain the final scaling law above, we 
have used that ^(k 7 ^ 0 ) ~ 77 “^ (k/ 0 ) 1 / y/N. We have neglected in the above scaling the possibility 
that fluctuations in rj may be large. 

Of course, the general framework can certainly account for this possibility, and one can directly plug 
in our ansatzes into (83) - however, we do not see general simplifcations that can be made, other than 
the crude scaling arguments here. Given our ansatzes, the expression in (131) is generally nonlocal and 
thus will not be as elegant as (126). Putting all of this together, we find (87). 

Essential to the scaling laws in (131) was that ^k~‘^ ~ above. However, this is not true in 
d < 2, where the sum will diverge at small k. We now give a heuristic argument that the typical scaling 
behavior we found above need not be parametrically different in these dimensions, consistent with what 
we found previously (in Section 3.2, the form of the conductivities is the same in all dimensions at weak 
disorder). To do so, we need to argue that there is a small modification that we can make to ip (so the 
contribution to the bound on conductivity from is qualitatively unchanged), yet which can remove 
the IR divergence from the viscous contribution. The natural guess is to modify (^(k) so that 7ij(k) = 0 
for |k|^ < (5, with <5 <C 1 a small constant. In this case, then we predict that 


1 ( N / 6 d = 1 

^ { Nlog(l/6) d = 2 


(132) 


and so if we choose a “reasonable” d (especially in d = 2 ), we can get acceptably finite viscous contributions 
to the conductivity. This can be accomplished by writing down 


(Q(k) + d-Q(k)) , 


(133) 


with (po given by ( 121 ), and q ^ Q, carrying no anomalous powers of N or 5. As typical elements of (p 
changed by an amount ~ 6, we expect that the conductivity in (126) has only changed by a small amount. 
Let us verify this is possible. We wish to find a solution to the highly overdetermined equation 

5ijQ(k)(l-6k,o)-Qo^Q(k)-^Q(k-q)^Q(q) = ^ Q(k-q)^(iQ(q), 0 < |k| < (134) 

The left hand side of the above equation scales as 1/y/N for typical disorder in Q. In the third term, 
this scaling follows from the lack of correlations between Q at different momenta. So long as we choose 
( 5 (k) = 0 in for |k|,^ < 25, then the “matrix elements” Q on the right hand side are also small. 
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We now look for Q(q) by solving a constrained optimization problem of the schematic form: 

a = Bx, (135) 

where a G and x G ni ~ dN <C n2 ~ (1 — 2'^6)N N, and B is a rectangular matrix, such 
that |x| is smallest, a is analogous to the (known) left hand side of (134), B is a known, truncated 
convolution-like matrix in a Fourier basis, and x plays the role of the undetermined, high momentum 
modes of (5(q) • S. This is a classic problem in constrained optimization with the solution [105] 

5Q(q) = a-(BBT)-^B. (136) 

Given that elements of a and B each scale as ^yl/N, we roughly estimate that <5Q(q) ~ 6, and so indeed 
it is possible to pick a small correction to ip which eliminates the divergence in the viscous contribution 
to the conductivity. Note that the matrix BB"*" is nearly diagonal (the off diagonal elements involve 
uncorrelated sums of random variables and so scale as l/\/iV instead of 1), and so there are no concerns 
about parametrically small eigenvalues of BB"*". 



Striped Models 


Thinking about resistivities turns out to be most convenient for models in d = 1, or with translational 
symmetry only broken in a single direction x, as noted in [11]. This follows from the general arguments 
that we make in Appendix C. Since we know that the current flow is a constant, to solve the variational 
problem exactly we need only vary (76) with respect to arbitrary v - the global minimum will corresponds 
to the true velocity field v. We hnd that v obeys the differential equation 


dx {rixxxxdxv) = [E {j^ - /u) 


(137) 


This second order linear differential equation cannot be solved exactly in general. We do not believe that 
closed form solutions exist in general for axx, though they can be found in special cases - for example, 
if is the only non-vanishing diffusive transport coefficient (as in [11], for a non-critical fluid), or if 
is the only non-vanishing coefficient, as in Section 4. In both of these cases, is not an invertible 
matrix - the zero eigenvector then provides a constraint which hxes Vx in terms of 

In particular, let us carry out this computation explicitly for the holographic striped models with 
equations of state given in Section 4. We must generalize the discussion to curved spaces, but this is not 
so difficult. The heat conservation equation implies (on a curved space) that 

j'^ = ^/yJ^^TSvx = e~^TSvx = constant. (138) 


Note that = 1, which simplifies calculations. After an appropriate generalization to curved space, we 
use (76) (on the true solution) to compute the inverse thermoelectric conductivity matrix: 


(ct“^)“^77“77^ = E 


v/7 {ixx {J^ - Qvxl^^Vif + 


(139) 


where 

(VfcUz + ViVk) - (140) 

For our set-up, parity symmetry in the y direction ensures that Vy = 0, and that the only non-vanishing 
components of are 

'JxxS'^'" = -lyyS^y = dxVx- (141) 
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Putting this together and using (61) to determine 77, S and Q, we obtain 





Sat 

4'^TSHtt 



2 


which gives (63). 


(142) 
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